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ABSTRACT 

The primary interest of the authors is in the area of grid 
generation, in paticular, optimal domain decomposion about realistic 
configurations. A grid generation procedure with optimal blocking 
strategies has been developed to generate multi-block grids for a 
circular- to -rectangular transition duct. The focus of this study is the 
domain decomposition which optimizes solution algorithm/block 
compatibility based on geometrical complexities as well as the physical 
characteristics of flow field The progress realized in this study is 
summarized in this paper. 


INTRODUCTION 

Most solution algorithms for irregular configurations presently solve 
a discrete form of the fluid equations of motion. The discretized 
equations must be solved on a discretized computational region, that is, 
on the grid nodes. A discretization process of grid generation is to 
establish a relationship between the physical and computational domain, 
thus allow the solution algorithm to be performed in the transformed 
domain. 

In principle it is possible to make a correspondence between any 
physical region and a single computational region. However, for general 
complicated three-dimensional configurations the resulting grid is likely 
to exhibit excessive skewness and coarseness. Despite the power and 
sophistication of present grid generation capabilities, it remains difficult 
to generate a reasonable, single-block grid about geometrically complex 
flow field configuration. 
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A better approach with complicated physical boundaries is to 
decompose the physical domains, each bounded by six curved sides and 
each of which transforms to a rectangular block in the computational 
domain(ref.l). This domain decomposition has the merits of flow solver 
efficiency, grid smoothness and orthogonality, and flexibility of gridding 
for complex geometries by choosing suitable block topologies. 

There are many approaches and philosophies which can be included 
under the domain decomposition strategies. Structured grid analysis of 
flow around complex geometries in widespread use is blocked 
decomposition, in which the solution domain js divided into regions 
with common internal boundaries. The grid lines at adjoining blocks 
can be set up to match in a variety of ways, with various levels of 
slope continuity. There is no theory which governs the way in which a 
flow field should be partitioned. A limited number of papers in the 
literature(ref.2-5) address basic flow-field-decomposition criteria and 
guidelines. Eiseman has developed an automatic block decomposition 
algorithm(ref.6). This algorithm is based on the geometrical 

complexities. However, no comprehensive, systematic studies have been 
done to determine the effect on the computed solution of using different 
blocks for the same geometry. Currently, the major bottleneck is in the 
design and implementation of the blocking plane(ref.7) 

To remedy this problem, the focus of this study is the domain 
decomposition which optimizes algorithm/block compatibility based on 
geometrical complexities as well as the physical characteristics of the 
flow field. 

Generally, the use of analytical shapes is not enough to satisfy 
the unusual geometrical requirement. Hence, sculptured curve/surface 
definitions such as Non-Uniform national B-Splines(NURBS), and quick 
elliptic grid refinement algorithms are developed. The application of 
these algorithms to grid adaption and domain decomposition is 
demonstrated. Grid generation associated with the circular-to-rect 
angular transition duct has been accomplished by applying these 
techniques. After careful consideration of the various alternatives, the 
structured muli-block approach must be chosen as the most suitable, 
from the point of view both grid generation and flow analysis. 
INS3D(ref.8), a three-dimensional multi-block incompressible Naver- 
Stokes code is used for this study. 
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DOMAIN DECOMPOSION PROCEDURE 


The work focuses on effects for optimizing domain decomposition 
strategies associated with circlar-to-rectangular transition duct The 
transition duct is designed to connect a typical circular engine exhaust 
to a high aspect ratio rectangular nozzle. The application presented is 
of considerable engineering importance in internal fluid flow designs. 

To take full advantage of the flexibility of multi-block structured 
techniques, one has to decide upon a suitable blocking topology to yield 
an optimal block arrangement for a given flow solver. 



Figure 1. Basic procedure for domain decomposition 


599 












This step is an art, requiring knowledge-based techniques and 
trial-and-error. Typically, a domain decomposition strategy is performed 
and then surfaces and volume grids are created for each individual 
block. Figure 1 shows the basic procedure. The desired characteristics 
of multi-block grid are strongly dependent on the needs of the flow 
solver and computer. 


GEOMETRY DEFINITION FOR TRANSITION DUCT 


The circle-to-rectangular transition duct can be designed by the 
equation of a superellipse. 
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Where a and b are the semi-major and semi-minor axes, respectively, 
of the superellipse. The x axis is coincident with the streamwise 
direction and the y and z axes are parallel to the major and minor axes 
of the duct. The cross-sectional area at a given streamwise location is 
given by 

r( — ) 2 

A a = 3 — <-f-)(4 ab) (2) 

r <T> 

Where ^ refers to the gamma fuction defined as 

r(t,)= r"( e ‘ r l )dt (n>0) (3) 

J 0 

The design procedure for the transition duct is to specify the 
streamwise variation of the semi- major and semi-minor axes, and 
superellipse exponent ti) defined by fifth order polynomial functions. 
Special cases of the superellipse include a circle( a=b, t)=oo) j an ellipse( 
a / b, ti=oo). 

Steady, incompressible, turbulent, swirl-free flow through a 
circle-to-rectangular transition duct has been studied experimentally 
(ref .9). For comparison, the same geometry has been simulated. Figure 
2 shows the lower half of the duct. The transition duct has an inlet 
diameter of 20.43cm, a length-to-diameter ratio of 1, and an exit plane 
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aspect ratio of three. The transition region length- to-diameter ratio is 
1.5. The cross-sectional area remains the same at the exit as at the 
inlet, but varies through the transition section to a maximum value 
approximately 15% above the inlet value. 


SURFACE REDISTRIBUTION AND REMAPPING 


The principal requirement for generating a grid about a complex 
geometry is the ability to efficiently redistribute points on the surface 
while maintaining the integrity of the geometry. Because the regions 
where high gradients are expected in the flow solution require high 
density of grid points. The redistributed surface grid is accomplished 
by evaluating the NURBS surface at the respective parameter 
associated with the desired distribution space. The NURBS is used for 
the standard surface description. The convex hull, local support, and 
variation disminishing properties of B- spline functions contribute to the 
generation of the well-distributed smooth grid. 

Let r=[xi(L 1 i), X 2 U. 11 ), xm^)} denote the parametric representation 
of the surface with coordinate (xi, xz, x 3 ) and parameter U, 11 ). A 
control point form of the NURBS surface(in 3-D) is defined as a tensor 
product formula in 4-D(ref.lO). 


B(i, = — - 

))H 0 

i=o;=o 


(4) 


Where B^(0 ; the M* degree B-spline basis functions in i direction. 
Bk*(Ti):the K* degree B-spline basis functions in j direction. 

Ha '■ weight (positive real values) 

Qa control points 

The control points of the surface are determined using tensor product 
formula associated with both f, and ^ parameters when a surface with 
a set of data point is given. It is called the inverse problem. 

The parametric space associated with NURBS is transformed as 
the normalized arc length distribution mesh. The original surface which 
is expressed into non-NURBS form can be converted to NURBS 
representation. It is important to note that the redistributed surface grid 
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is obtained by evaluating the NURBS surface at the desired distribution 
point. However, the NURBS surface is smooth and has kept precisely 
the fidelity of the original surface. A redistributed surface resulting 
from the NURBS surface is demonstrated in Figure 3. 


ELLIPTIC REFINEMENT 

The grid should exhibit the desirable qualities of smoothness and 
orthogonality especially near the solid surface. Some results of algebraic 
grid methods such as transfinite interpolation functions can not meet 
the requirement of the grid qualities on certain geometries. Thus, 
elliptic grid methods are used to make up the shortcoming of the 
algebraic methods(ref.ll). 

Derivation of the control function based on grid metrices begins 
with analyzing the elliptic grid equation. Define the physical and 
computational space with r=(x i, xz, *3) and 12= U 1 , £ 2 , 4 3 ), respectively. 
The covariant and contravariant vectors then appear as follows: 
at = co variant base vectors z? (« = 1, 2, 3) 

a = contravariant base vectors (i=l, 2, 3) 

gu = a • a. i = g* U=l, 2, 3), 0=1, 2, 3) 

g u = d • d = g * (i=l, 2, 3), 0=1, 2, 3) 

g = det\gij\ = [ ai • ( U 2 X m)] z 
( gu ) = derivative of g a with respective to t, k 

= m (i = l, 2, 3), 0=1. 2, 3), (fc= 1, 2, 3) 

The elliptic grid generation system used in this study (ref. 12) is 

+ Tg^PkZ? = 0 (5) 

t = I/=I k = i 

The determinations of the three control functions for the general 
three-dimensional case can be summarized as follows. The three 
components of the elliptic grid generation equation(5) provide a set of 
three equations, 

( x t ‘), Pk = -jtjky ( zw>» /=i. % 3 (6) 
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that can be solved simultaneously at each point for the three control 
functions, Pk, k= 1 , 2 , 3. The derivatives here are represented by 
central differences. The smoothness is established by replacing the 
control function at each point with the average of the four neighbors in 
the two curvilinear directions other than that of the function. Thus, 

pm 1 , t 2 , £ 3 )=-^-[PiU\ t 2 +i, s 3 ) + Pitt 1 , e-i, t 3 ) 

+ Pitt 1 , i, 2 , ^ 3 + l) + Pitt 1 , (, 2 , (7) 

with analogous equations for P 2 and P 3 . No smoothing is done in the 
direction of the function because to do so would smooth the spacing 
distribution. The use of smoothed control functions evaluated from the 
algebraic grid produces a smooth grid that retains essentially the 
spacing of the algebraic grid(ref,13). 

An application of these control functions results in a smooth-nearly 
orthogonal grid in fewer iterations of the elliptic solver. These control 
functions are applied in surface/volume grid refinement. It can be 
observed that the elliptic grid provides smoothness and near 
orthogonality in Figure 4. 


BLOCK INTERFACES 


There are three distinct configurations along the axial direction 
which are a constant diameter circular entrance section, the transition 
section, and a rectangular extension section. Because of the drastic 
changes in the flow direction, a multi-block grid topology has been 
adapted. This provides a smooth discretization of the entire volume 
inside the transition duct 

Each block has its own curvilinear coordinate system irrespective 
of that in the adjacent blocks. In order to keep complete continuity of 
grid lines across the interfaces between blocks, the linkage among the 
various blocks can be set by interpolation or fixing an actual boundary. 
The grid points on an interface of one block are coincident in physical 
space with those on another interface of the same or another block. 
This facilitates the interface of block treatment without an application 
of interpolation. This philosophy allows the flow solver to be run on 
the multi-block grid system. 
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APPLICATION 


The two previously described methods, elliptic and algebraic, are 
merged appropriately to obtain various multi- block grids on the 
transition duct configuration. Moreover, by utilizating an approprite 
blending of these methods allows a grid generation capability ranging 
from cases with only a few blocks to cases with hundreds of blocks 
for various shapes. 

To demonstrate the capability of the present procedure, several 
types of multi- block grids are designed such that the grid cell aspect 
ratio and orthogonality are maintained with a reasonable range for all 
the duct cross section. Figure 5-8 present several examples. For 
comparative purposes, all types of domain decomposition are 
implemented in the same flow solver code using equivalent boundary 
conditions. Even the flow conditions for the computational analysis are 
chosen to correspond to the experiments reported in Ref. 8. Inlet 
conditions are as follows: 

Reynolds number = 390,000(based on U and D) 

Bulk velocity(U) = 29.95 m/s 
Core velocity(Uc/U) = 1.083 
Friction velocity(Utau/U) = 0.04063 


RESULTS AND DISCUSSION 

The results of the flow calculation are shown in Figures 9-12. 
Pressure contours of the axial flow component are plotted in Figure 9 
(a)-(b) which are on the x-y plane and x-z plane, respectively. The 
velocity contours on both planes are demonstrated in Figure 10.(a)-(b). 
The transition segment produces saddle-shaped pressure distributions in 
the y-z plane as shown in Figure 11. Within the boundary layer the 
velocity is reduced as shown in Figure 12, but the cross-stream 
pressure gradient imposed by the flow outside the boundary layer is 
not reduced. This can result in significant flow turning in the boundary 
layer and is refered to as skew-induced secondary flow or cross flow. 
The rate of cross flow production increases as the amount of 
streamline curvature increases(ref,13). Figure 13 shows the comparison 
of sidewall normalized static pressure between the experimental data 
and the computed solutions in the upper quarter of the duct only. The 
results of computational simulation and the measured data of Ref. 9., 
show a very good agreement. 
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CONCLUSIONS 


Several blocking strategies have been considered in order to 
analyze the transition duct flow. Optimal domain decomposition is 
dependent on the method of flow solver and flow properties as well as 
geometry concerns. This domain decomposition gives us the following 
merits: 

• A higher solution accuracy and faster convergence for the 
computational fluid dynamics solver. 

• Savings in the CPU time and the amount of memory. 

• Maximization of the grid quality and optimization of the grid 
distribution. 
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STATION 1 


Figure 2. The lower half of the duct 
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Figure 3. A NURBS surface 



Figure 4. Elliptic grid refinement 
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Figure 5. Two-block system along the flow direction 



Figure 6. Two blocks, O-type grid system. 
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Figure 7. Two blocks, O-H type grid system. 



Figure 8. 5 blocks FI -type grid system. 
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Figure 9. Pressure contours (a)x-y plane, (b)x z plane. 


(a) 



(b) 



Figure 10. Velocity contours (a)x-y plane, (b)x-z plane. 
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Figure 11. Vorticity contuors. 



Figure 12. Boundary layer effects. 
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peri peri al wall along the upper quarter. 


Figure 13. Comparison of peripheral wall pressure 
coefficient distribution between computed and measured 
data, (symbol: experiment, simbol- line '-calculation) 
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Figure 13. Comparison of peripheral wall pressure 
coefficient distribution between computed and measured 
data (symbol: experiment, simbol-line:calctilation) 


613 




NEW METHODS/APPROACHES/ 
APPLICATIONS (1) 
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